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ABSTRACT 

'St 

Four  modes  of  possible  rotor  motion  are  defined:  synchronous  preces¬ 
sion  (without  rub),  partial  rub,  full  annular  synchronous  rub  and  reverse 
whirl.  A  model  is  developed  to  investigate  the  relationships  between  a 
rotational  system's  characteristic  parameters  and  the  conditions  of  equili¬ 
brium  for  the  full  annular  rub  motion.  The  dynamics  of  an  unconstrained 
rotor  with  unbalance  are  first  established.  The  impact  of  constraining  the 
rotor's  motion  by  a  rigid  casing  with  a  finite  clearance  is  then  analyzed. 
For  any  given  system,  the  speed  regimes  in  which  full  annular  synchronous 
rub  is  in  equilibrium  are  defined  as  a  function  of  the  system's  unbalance, 
clearance,  rotor  stiffness,  damping  and  the  coefficient  of  friction  between 
the  rotor  and  casing.  The  rigid  stator  restriction  is  subsequently  allevi¬ 
ated  and  the  equilibrium  conditions  for  full  annular  synchronous  rub  are 
reanalyzed.  Additionally,  reverse  whirl  as  a  possible  stable  mode  of  rotor 
motion  is  studied.  Specific  rotational  systems  are  then  analyzed  according 
to  the  model  developed  to  demonstrate  tne  effect  of  varying  particular 
parameters.  .  ■} 
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INTRODUCTION 


All  rotating  machinery  have  some  practical  degree  of  unbalance  which 
drives  flexural  vibrations  in  the  system.  These  vibrations  may  be  readily 
calculated  in  terms  of  the  system's  mass,  unbalance,  rotational  speed, 
radial  stiffness  and  damping  imposed  by  the  medium  in  which  the  rotating 
element  operates.  Should  the  magnitude  of  the  vibrations  exceed  the 
physical  constraints  of  the  enclosure  which  houses  the  rotating  element,  the 
potentially  destructive  phenomena  described  here  as  rotor  rub  occurs. 

Shaft  vibrations,  misalignment,  blade  creep  and  thermal  expansion  are 
some  of  the  physical  realities  which  precipitate  design  of  a  radial 
clearance  between  the  high  speed  rotor  and  surrounding  casing  in 
turbomachinery.  Theoretical  efficiencies  of  the  turbomachinery  are  reduced 
though,  by  blade  tip  losses  which  result  from  the  presence  of  the  radial 
clearance  between  the  rotor  blade  tip  and  the  casing.  Demand  for  improved 
system  efficiency  thus  drives  a  reduction  of  the  blade  tip  clearance  and  so 
increases  the  likelihood  of  rotor  rub.  Rotor  rub  in  fact  may  be  a  design 
condition  and  so  accounted  in  the  design  by  the  use  of  abradable  materials. 
[1] 

Prior  to  the  occurrence  of  rub,  the  rotating  element’s  dynamics  may  be 
described  as  synchronous  precession,  driven  by  the  system's  residual 
unbalance.  If  rub  develops,  it  may  occur  over  a  fraction  of  the  rotor's 
precessional  orbit  (termed  partial  rub),  or  may  persist  over  the  full  orbit 
(full  annular  rub).  Of  these  two  conditions,  full  annular  rub  is 
potentially  more  destructive  and  is  the  subject  of  this  thesis.  For  further 
information  regarding  partial  rub,  attention  is  drawn  to  reference  [2-A]. 
Examples  of  full  annular  rotor  rub  occur  with  rotor  blade-to-casing  or 


rotor-seal  rub,  as  well  as  with  a  rotating  shaft  having  clearance  in  a 
bearing.  [5]  The  rub  introduces  a  reaction  force  at  the  rotor's  periphery 
which  strongly  influences  the  rotor's  dynamics  and  may  additionally  excite 
the  casing  (further  referred  to  as  the  stator)  into  motion.  The  reaction 
force  may  be  broken  down  into  its  normal  and  tangential  (friction)  force 
components;  the  two  being  linearly  related  by  the  coefficient  of  friction. 
Friction  may  be  moderated  by  proper  lubrication  (such  as  in  the  case  of 
bearing  rub)  as  well  as  through  suitable  selection  of  materials.  Other 
design  considerations  or  abnormal  operating  conditions  may  preclude  these 
preventive  measures  though,  and  so  friction  may  take  on  a  particularly 
destructive  nature. 

Friction  is  the  resistive  force  which  acts  tangentially  at  the 
interf-.ce  Detween  two  surfaces  in  relative  motion.  By  acting  tangentially 
o:.  t:  -  rot  Jr's  periphery,  it  essentially  establisnes  a  torque  which  counters 
tr :•  :r:  ving  input)  torque  for  the  rotating  element.  The  result  is  either  a 

-  :  .  r ,  t  j'  .  r , . >  1  sped  or  d  higher  input  torque  required  to  maintain  a 

a*  The  effect  of  friction  on  the  mating  surfaces  is  to  cause  a 

.  , .  ; .  o-jt  bu.ioup  and,  in  the  very  least,  some  degree  of  grinding  wear  and 
deformation  of  the  rotor  and  seal  or  bearing  annulus.  The  more  extreme  case 
encompasses  catastrophic  destruction  of  the  machine.  [2]  Typically 
accompany ir.g  these  mechanisms  is  the  annoying  and  often  unacceptable 
acoustic  phenomenon  due  to  frictional  vibration  instabilities  known  as  stick 
slip.  Dynamically,  the  occurrence  of  rotor  rub  may  destabilize  the 
synchronous  precession  of  the  rotor.  If  friction  forces  are  of  large  enough 
magnitude,  the  counter-torque  developed  may  instigate  a  reverse  whirl 
condition  wherein  the  precession  is  in  a  direction  opposite  to  the  direction 
of  rotation  of  the  shaft.  The  nature  of  reverse  whirl  is  such  that  the 


LO- 


1 .  DEVELOPMENT  OF  THE  MODEL 

The  characteristics  of  rotating  machinery  designed  for  the  generation 
or  transmission  of  power  vary  greatly  with  their  particular  application.  An 
exact  analysis  of  the  motion  is  extremely  complex,  but  many  meaningful 
results  can  be  determined  through  employment  of  reasonable  simplifying 
assumptions  in  a  single  degree  of  freedom  model.  Even  with  these  simplifi¬ 
cations  though,  numerous  variables  are  necessary  to  describe  the  motion 
involved. 

The  model  in  this  thesis  is  displayed  in  Figure  1.  The  rotating 

element  consists  of  a  central  disk  of  mass,  M  ,  mounted  symmetrically  on  a 

laterally  flexible  but  torsionally  rigid,  massless  shaft  of  stiffness  K  . 

The  shaft  stiffness  is  determinable  knowing  the  shaft's  material  type, 

cross-sectional  properties  and  length.  Tne  shaft  is  then  supported  at  each 

end  by  rigid  bearings  without  clearance.  The  natural  frequency  of  the 

rotor,  r >  the  frequency  of  excitation  at  which  resonance  occurs,  is  given 

by  =  /K  / M^ .  The  rotor  disk  is  enclosed  in  its  casing,  separated 

radially  by  a  symmetric  clearance,  r  ,  existing  when  the  rotor  is  at  rest. 

(The  effects  of  gravity  are  neglected  in  this  analysis).  The  casing,  or 

stator,  is  character ized  by  its  mass,  M  ,  and  its  mounting  stiffness,  K  . 

s  s 

As  in  the  case  of  the  rotor,  the  natural  frequency  of  the  stator  is  given  by 
-  /iWM^.  Damping  is  introduced  for  motion  of  the  rotor  and  stator 
respectively  by  the  damping  coefficients,  D^  and  D  .  It  is  the  presence  of 
damping  in  both  the  rotor  and  stator  which  limits  the  magnitude  of  their 
responses  at  their  respective  resonant  frequencies.  Damping  also  causes  a 
phase  lag  between  the  system's  excitation  and  response.  For  ease  of 


analysis,  the  model  has  been  "compressed"  in  Figure  2  to  a  planar  diagram 


with  the  element's  dynamic  characteristics  represented  as  in  a  simple  linear 
resonator  system:  a  spring-mounted  mass  with  a  dashpot  in  which  the  energy 
of  the  motion  is  dissipated.  Both  radial  and  longitudinal  symmetry  are 
implicit  in  this  model. 


2.  FREE  ROTOR  DYNAMICS 


It  is  assumed  that  under  design  conditions,  the  rotor  is  driven  at  a 
constant  speed,  uj,  independent  of  the  external  load  applied.  This  assump¬ 
tion  essentially  neutralizes  the  effects  of  fluctuations  in  the  design  load 
due  to  damping  or  rub  related  friction.  Under  ideal  conditions,  it  is 
imaginable  that  the  rotor  is  perfectly  balanced.  Thus  when  driven  at  its 
design  speed,  the  rotor  center  remains  at  the  center  of  rotation;  aligned 
with  the  axis  defined  by  the  supporting  bearings'  centers.  This  idealiza¬ 
tion  is  not  realized  in  practice  though,  and  to  some  extent,  a  residual 
unbalance  exists  in  the  mass  distribution  of  the  rotor.  This  unbalance  may 
be  depicted  in  the  planar  model  by  locating  the  mass  center  (or  center  of 
gravity),  g,  at  a  distance  6  from  the  center  of  the  disk.  The  effect  of  the 
unbalance  is  to  establish  a  centrifugal  force  acting  upon  the  center  of 
gravity  while  the  rotor  is  in  motion.  This  centrifugal  force,  which  is 
proportional  to  the  square  of  the  rotational  speed,  acts  as  the  driving 
force  for  the  rotor's  dynamical  motion. 

An  analysis  of  the  motion  of  the  free  rotor  ( i . e . ,  without  displacement 
constraints  imposed  by  the  stator)  provides  the  baseline  condition  for  which 
the  rub  condition  may  be  compared.  The  equation  of  motion  for  the  free 
rotor  is  derived  by  applying  Newton's  2nd  Law  to  the  rotor: 

(sum  of  forces  acting  on  the  rotor)  = 

(mass  of  the  rotor)  x  (acceleration  of  the  rotor  center  of  gravity). 
Referring  to  the  free  body  diagram  for  the  rotor,  Figure  3>  the  forces  which 
must  balance  in  the  steady  motion  of  the  rotor  are  described  in  the 


following : 
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a.  Centrifugal  Force  -  As  described  above,  this  force  results  from  the 
presence  of  a  residual  unbalance;  i.e.,  the  center  of  gravity  does  not 
coincide  with  the  center  of  rotation.  The  centrifugal  force  acts  on  the 
center  of  gravity  and  is  directed  radially  outward  from  the  bearing  axis 
(stationary  reference  position  about  which  the  rotor  precesses). 

b.  Spring  Force  -  The  spring  (or  restoring)  force  acting  radially 
inward  on  the  rotor  center  results  when  the  rotor  center  is  displaced  from 
the  bearing  axis.  It  is  assumed  in  this  model  that  the  spring  force  is  a 
linear  function  of  the  shaft  stiffness  (assumed  constant)  and  magnitude  of 
displacement. 

c.  Damping  Force  -  The  damping  force  acts  on  the  rotor  center  such 
that  it  resists  translational  motion  by  the  rotor.  The  damping  force  is 
assumed  to  be  a  function  of  the  damping  coefficient  experienced  by  the  rotor 
in  its  particular  medium  and  the  linear  component  of  the  rotor's  velocity. 

The  nature  of  the  rotor's  dynamics  readily  lend  the  analysis  to 
description  using  complex  variables.  Defining  zq  as  the  time  dependent 

4 

position  of  the  rotor's  center,  referenced  such  that  z^  =  0  coincides  with 
the  rest  condition  where  the  rotor  center  is  aligned  with  the  bearing  axis, 
then  the  equation  of  motion  may  be  written  as: 

IF  =  [-D  z  -  K  z  ]  -  M  z  (1) 

r  c  r  c  r  g 


where  z  (t)  =  the  vector  position  of  the  center  of  gravity  of  the  rotor. 
S 

Furthermore,  z  (t)  may  be  broken  down  into  the  vector  sum: 

8 


z  (t)  =  z  (t)  +  z  (t) 
g  c  gc 


where  z  (t)  represents  the  vector  position  of  the  rotor  center  of  gravity 
gc 

relative  to  the  disk  center. 
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The  equation  of  motion  then  becomes: 
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M  (z  +z  )=-Dz  -Kz 
r  c  gc  r  c  r  c 


or, 


Mz  +  Dz  +  Kz  =*  M  z 
rc  rc  rc  r  gc 


(2) 


z  (t)  is  simply  a  function  of  the  radius  of  unbalance  and  rotational 
speed: 


z  (t)  »  6e'^u>t 
gc 

Assuming  synchronous  whirl  response,  the  motion  of  the  rotor  may  be 
described  by: 


Making  these  substitutions  into  Equation  (2),  and  simplifying,  results  in 
the  following: 


+  LM  6]u>2ejult 

2  (t)  - - 2 - 

[(Kr  -  u  Mp)  ♦  J(uDp)] 


This  harmonic  motion  of  the  rotor  is  termed  forward  synchronous  precession 
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since  the  orbital  motion  of  the  rotor  is  characterized  by  the  forcing 
function's  angular  frequency  and  direction.  Breaking  down  the  rotor's 
mechanical  impedance, 


C(Kr  -  u>  Hp)  ♦  J(uDp)]  , 

into  its  magnitude  and  phase  then  leads  to  the  harmonic  response  equation: 


[M  6]a)2e^(wt-<^) 

zcU)  = - -y2 

[(K  -  uj2M  ) 2  ♦  ( uD  ) 2] 

r  r  r 


(3) 


The  phase  angle  <p  describes  the  lag  between  the  rotor's  response  and  the 


forcing  function.  This  phase  lag  which  is  due  to  damping  is  given  by: 
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It  is  convenient  to  introduce  non-dimensional  variables  into  Equation 
(3)  for  ease  of  future  analysis.  By  defining  a  non-dimensional  rotor  speed 
and  damping  ratio  as  given  below: 

0  =  w/wNr 


?  = 


2/K  M 
r  r 


then  Equation  (3)  may  be  rewritten: 

■»  _  _602  _  ^j(u)t-<t>) 

2°  =  "  2  2  2  1/2 
t (1  -  B  )  ♦  (2c6)  ] 

Similarly,  the  phase  relationship  may  be  rewritten: 

r1  -U562- 


(JO 


4)  =  tan 


(4a) 


(1  -  Bc) 


The  plots  of  zc  vs.  0  and  <j>  vs.  0  are  typical  amplitude  and  phase  plots 
for  a  second  order  system,  such  as  those  in  Figure  4. 

Consider  the  case  of  light  damping,  r,  <<  1.0.  At  low  speeds,  i.e., 

.■+  .  2 

0  <<  1.0,  the  magnitude  of  Equation  (4)  reduces  to  |zc|  *  60  .  As  the  speed 

increases,  the  displacement  increases  non-linearly  until  the  maximum  dis- 

2  1/2 

placement  occurs  at  0  -  (1  -  2^)  : 


c '  max 


1  /  2 


2C< 1  -  C  > 

Here  it  can  be  seen  that  the  maximum  displacement  is  solely  a  function  of 


.V.'  .* 
_  V 


/Sjw 

m 


.. 

»: 

•.v.y; 

..-.V 


P 

'•'VI 


4  s 


\'s: 


a  .* 


r 


\.\.N 

•v-.n 

.'I 

V  .’J 

agasjjf 


the  radius  of  unbalance  (directly)  and  the  damping  ratio  (inversely). 

Further  increases  in  speed  bring  a  reduction  in  displacement.  For  very  high 
speeds,  i.e.,  6  >>  1.0,  the  displacement  approaches  6  (a  constant). 


The  phase  relationship  for  light  damping  is  given  by  Equation  (4a). 

For  low  speeds,  8  <<  1  .0,  <f>  is  approximately  zero.  Essentially,  the  forcing 
function  and  the  response  are  exactly  in  phase.  As  the  speed  increases 
though,  and  approaches  the  resonant  speed,  the  damping  force  becomes 
significant  and  the  response  begins  to  lag  the  forcing  function.  Near 
resonance,  a  rapid  phase  shift  occurs  such  that  at  resonance,  <J>  =  90°.  As 
the  speed  further  increases,  4>  asymptotically  approaches  tr.  In  the  limit, 


i.e.,  for  t;  =  0;  $  =  0  for  all  values  of  8  <  1.0  and  $  =  it  for  all  values 


of  8  >  1.0. 


3.  CONSTRAINED  ROTOR  DYNAMICS 


As  previously  described,  rotating  machinery  designed  for  the  generation 
or  transfer  of  pow^r  are  typically  enclosed  within  a  casing  which  minimizes 
the  clearance  between  the  rotating  element  and  casing.  The  dynamics  of  the 
constrained  rotor  motion  may  simply  reduce  to  those  of  the  free  rotor  as 
presented  in  the  preceding  section  if  the  rotor  and  stator  never  come  into 
contact.  Tnis  no-rub  condition  exists  if  either  of  two  basic  criteria  are 
in  effect: 

(1)  Tne  maximum  amplitude  of  free  rotor  displacement,  |z  I  ,  is  less 

c  max 

than  the  constrained  rotor  clearance,  or 

(2)  The  rotor  is  operated  within  a  speed  regime,  characterized  by  a 
speed  w,  such  that  tne  amplitude  of  free  rotor  displacement,  ] z  ( w) | ,  is 
less  than  the  constrained  rotor  clearance. 

As  will  become  evident  through  the  following  analysis,  condition  (2)  above 
should  be  further  restricted  to  speed  regimes  where  u>  is  less  than  the  rotor 
natural  frequency,  These  two  cases  are  depicted  in  Figure  5  using  the 

plot  of  amplitude  vs.  speed. 

If  in  fact  the  rotor  and  stator  come  into  contact,  the  dynamics  may 
become  quite  complex.  The  analysis  will  be  first  conducted  with  the 
assumption  that  the  stator  is  perfectly  rigid  (though  this  restriction  will 
be  alleviated  later)  and  only  full  annular  rub  will  be  considered. 

The  free  body  diagram  for  the  rub  case  is  presented  in  Figure  6.  Prior 
to  contact,  Equation  (4)  describes  the  rotor's  motion.  It  can  be  seen  that 
upon  the  initiation  of  rub,  the  rigid  stator  limits  the  rotor's  displacement 
while  introducing  a  reaction  force,  R(t),  at  the  point  of  contact.  This 


displacement  rotor  displacement 
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reaction  force  is  comprised  of  a  normal  force  and  a  tangential  (friction) 
force.  Since  the  displacement  is  constant,  the  restoring  force  acting  on 
the  rotor  is  also  constant  (i.e.,  independent  of  rotor  speed,  w) .  As  in  the 
case  of  the  free  rotor,  the  equation  of  motion  for  the  rub  condition  may  be 
developed  from  Newton's  Second  Law: 


M  z  =~Dz  -Kz  -R(t) 
r  g  r  c  re 


(6; 


Jwt 


Substituting  z  =  (z  +  z  ),  and  R(t)  =  (N  +  ipN)e  ,  then  Equation 
g  gc  c 


(6)  becomes  (after  rearranging  terms): 


Mz  +  D  z  +  K  z 
re  r  c  r  c 


-M  z 
r  gc 


R(t) 


Jut 


Since  z  =  6e  ,  then  the  above  equation  may  be  solved  for  z  by 
gc  c 


assuming  synchronous  precession: 


„  [M  6u>  -  (N  +  juN)] 

2  (t)  = - - eJu 

[(K  -  uTM  )  +  j(a>D  J] 

r  r  r 


(7) 


By  dividing  through  Equation  (7)  with  the  radial  clearance,  r  ,  and 


utilizing  the  previously  defined  non-dimensional  speed  and  damping  ratio 
variables,  Equation  (7)  is  transformed  to  its  non-dimensional  form: 


(t) 


_Lc6__3_^n_  +  _jijn2J_  ejoot 

[(1  ~  S2)  +  JC2C8)] 


(3) 


where  the  new  variables  are  defined: 
e  =  6/r 


n  =  N  /  k  r 


The  right  hand  side  of  Equation  (8)  is  reduced  to  a  magnitude  and  a  phase  in 
the  following  rotor  rub  equation: 

1/2 


z  ( t) 
c 


£  ( e62g_n)2_i_(jjn)2] 


J(wt  -  i>  -  y) 


td  -  s2)2  +  (2;e)2] 


1/2 


(9) 
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where 


.  -1  (yn) 

Y  =  tan  — ~%~L — 

(eg  -  n) 

The  following  observations  may  be  drawn  regarding  the  physical 
implications  of  the  rotor  rub  equation.  The  parameter  6  is  a  measure  of  the 
rotor's  unbalance,  and  as  such  is  not  designed  into  the  rotor,  but  results 
from  manufacturing  processes  which  are  unable  to  meet  with  ideal  design 
conditions.  The  value  of  though  is  designed  into  the  system,  and  so 
tolerances  must  be  placed  on  5  in  order  to  ultimately  control  the  value  of 
e,  a  measure  of  the  rotor's  displacement-to-clearance  ratio.  In  the  case  of 
the  free  rotor,  any  increase  in  6  leads  to  a  direct  increase  in  the  rotor's 
maximum  displacement.  Since  the  maximum  displacement  is  now  constrained 
though,  a  normal  force  develops  at  the  interface  between  the  rotor  and 
stator  when  the  rotor's  runout  equals  the  clearance.  The  maximum  normal 
force  developed  thereafter  is  directly  related  to  the  value  of  e.  Referring 
to  Figure  5,  the  normal  force  developed  for  rotor  speeds  less  than  the  rotor 
resonant  speed  may  be  crudely  estimated  as  the  magnitude  of  force  necessary 
to  compress  a  'spring'  of  stiffness  equivalent  to  the  rotor's  impedance  a 
distance  equal  to  tne  difference  between  the  free  rotor's  displacement  and 
the  radial  clearance.  For  speeds  greater  than  the  rotor  resonant  speed 
(6  >  1.0)  the  characteristics  of  the  combined  rotor-stator  govern  the 
system's  behavior  and  a  more  detailed  analysis  is  necessary. 

Meanwhile,  a  frictional  force,  f  =  yn,  develops  at  the  rotor-stator 
interface  in  opposition  to  the  relative  motion  between  the  rotor  and  stator. 
The  coefficient  of  friction,  y,  depends  primarily  on  the  material  types  in 
contact  and  the  degree  of  lubrication  present.  For  unlubricated  metal-to- 
metal  applications,  typical  values  of  y  may  range  from  0.1  -  0.3.  [7]  The 
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frictional  force  may  be  replaced  by  a  force  and  a  couple  acting  upon  the 

rotor's  center.  The  torque  thus  created  acts  as  an  additional  load  upon  the 

rotor's  driving  force.  Equation  (9)  emphasizes  the  phase  lag,  (41  +  Y), 

between  the  response  of  the  rotor  and  the  forcing  function.  4  was  defined 

previously  as  the  phase  lag  caused  by  the  damping  force  in  opposition  to  the 

linear  component  of  the  rotor's  velocity  and  was  given  by  the  expression, 

-1  2 

<(>  =  tan  (2cB)/(1  _  8  )•  Y  is  the  phase  lag  caused  by  the  frictional  force 
which  further  opposes  the  rotor  motion  as  previously  given  by: 
y  =  tan  1 (yn)/( e82  -  n) . 

Since  rub  is  presumed,  the  magnitude  of  the  rotor's  displacement  is 
exactly  equal  to  the  radial  clearance  (given  a  rigid  stator).  This 
condition  results  in  the  magnitude  of  Equation  (9)  being  equal  to  unity. 

The  independent  variable  in  the  equation  is  the  rotor  speed.  The  values  of 
y ,  z,  and  e  are  assumed  constant,  and  so  the  values  of  the  normal  force  and 
phase  angle  may  be  solved  for  any  value  of  B.  This  arises  from  the 


condition: 


r  ,  02  .2  ,  .2-.1  /2 
Lis5_-_n)  t  lt!n!  ]  ,  n 

1/2  ,U 

[(1  -  B2)2  +  (2tB)2] 


.-.y.v 

V>’v 

/%'vJ 


Following  a  few  algebraic  steps,  the  following  quadratic  equation  in  n  is 


derived: 

(1  +  y2)n2  +  2(1  -  B2  +  2;By)n  +  [(1  -  82)2  +  (2c6)2  -  e2 B**]  =  0  (10) 


which  when  solved  yields  the  relationship  for  n  as  a  function  of  B: 


(1  +  y2) 


Equation  (11)  indicates  two  possible  solutions  for  the  normal  force. 
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It  is  the  upper,  'positive'  root  which  is  of  primary  interest,  yielding  a 
condition  of  n  «=  0  which  is  consistent  with  the  rotor  speed,  B,  at  which  rub 
initiates.  As  B  varies,  so  too  the  values  of  the  normal  force  and  phase 


angle  vary  together  as  necessary  to  maintain  a  balance  of  all  forces  acting 
on  the  rotor. 

It  is  evident  from  Equation  (11)  that  throughout  the  regime  of  possible 
rotor  speeds,  depending  upon  the  values  of  the  system  parameters  y,  £  and  e; 
the  normal  force  may  be  positive,  negative  or  imaginary.  A  negative  value 
of  the  normal  force  here  implies  tensile  stresses  existing  between  the  rotor 
and  stator.  This  is  essentially  infeasible,  and  so  for  all  values  of  n  <  0, 
the  rotor  and  stator  are  deemed  not  to  be  in  contact.  The  speed  at  which 


rub  is  actually  initiated  is  determinable  by  solving  Equation  (11)  for  the 
value  of  0  at  which  n  =  0: 


(e2  -  1) 


This  value  of  B  then  results  in  |zc|  equal  to  rc  in  both  Equations  (A)  and 
(9).  It  should  be  immediately  noted  that  if  e/2 c  <  1.0  for  a  particular 
system,  then  rub  in  fact  does  not  occur  at  any  rotor  speed.  The  inequality 
above  is  simply  the  numerical  equivalent  to  a  previously  described  condition 
-  i.e.,  the  maximum  amplitude  of  the  free  rotor  displacement  is  less  than 
the  radial  clearance  allowed  the  rotor.  This  is  indeed  the  most  common 
de3igned-for  condition.  Through  service  use  though,  the  value  of  e  for  a 
rotating  machine  may  expectedly  change  since  both  the  degree  of  unbalance 
and  the  radial  clearance  are  affected  by  such  physical  realities  as  creep 
caused  by  high  centrifugal  loading,  thermal  and  aerodynamic  stresses  caused 
by  interaction  between  the  rotor  and  the  working  fluid,  or  possible  shaft 
bowing  due  to  uneven  heating  or  cooling  or  simple  static  loading. 


It  is  the  intent  of  this  analysis  to  determine  the  behavior  of  rotating 
machinery  as  the  rotor  speed  is  increased  beyond  the  rub  speed  for  systems 
where  e/2 s  >  1.0.  Figure  5(b)  provides  a  simple  plot  of  the  amplitude  vs. 
speed  relationship  for  the  rotor.  It  is  known  that  for  6  <  8  Equation 
(M)  describes  the  rotor's  motion.  Furthermore,  for  rotor  speeds  within  the 
regime  Bpub  <  8  <  0' ,  it  is  assured  that  the  rotor  and  stator  will  be  in 
contact  with  the  normal  force  given  by  Equation  (11).  For  rotor  speeds 
greater  than  6'  though,  Figure  5(b)  can  only  provide  basic  insight  into  the 
problem.  It  was  alluded  earlier  that  while  the  rotor  is  undergoing  a  full 
annular  synchronous  rub,  then  it  is  the  characteristics  of  the  combined 
rotor-stator  which  govern  the  system's  performance.  This  performance  may 
best  be  analyzed  by  investigating  the  normal  force  and  phase  angle 
experienced  by  the  rotor. 

Equation  (11)  for  the  normal  force  may  be  rearranged  algebraically  to 
yield  the  following  expression: 


-  c(i-62)  ♦  u(2Cb>]  ±  (i+p2)1/2/£62)2  -  Mizeh-i.iui2 

n  - - z - U-t-U-l—  (12) 

(1  ♦  U  ) 

The  three  possible  regimes  of  n  (positive,  negative  or  imaginary)  may  be 
examined  through  this  equation.  It  is  instrumental  to  first  analyze  F,  the 
sum  of  the  external  forces  acting  on  the  rotor  center,  by  resolving  it  into 
orthogonal  components  (see  Figure  7)  _  those  forces  acting  parallel  to  zc 
(designated  the  radial  or  r-direction)  and  those  forces  acting  normal  to 
zc  (designated  the  tangential  or  t-direction) . 

4 

F  is  further  resolved  into  a  component  due  to  rotor  displacement  and 
a  component  due  to  the  mass  eccentricity  itself  in  Figure  8. 

By  combining  the  non-dimensional  relationships  in  Figures  (7)  and  (8), 
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a  geometric  interpretation  for  the  relationship  between  forces  in  Equation 
(11)  is  attained  in  Figure  9. 

By  applying  the  Pythagorean  Theorem  to  Figure  9,  the  variables  become 
related  by: 

( eS2)2  »  (2cB  +  un)2  +  [(1  -  B2)  +  n]2  (13) 

Physically,  this  states  that: 

[the  component  of  centrifugal  force  acting  on  the  rotor  center  due  to  mass 
2 

eccentricity],  e3  ,  is  balanced  tangentially  by  [the  damping  and  friction 
forces],  and  radially  by  [the  normal  and  elastic  forces  lessened  by  the 
centrifugal  force  acting  on  the  rotor  center  due  to  rotor  displacement].  By 
solving  for  the  normal  force  itself,  Equation  (13)  yields: 

1/2 

n  -  -(1  -  B2)  ♦  { ( cfi2)2  ~  (2;B  +  un)2}  (14) 

The  parallel  between  this  equation  and  Equation  (12)  is  immediately 
apparent.  The  following  generalizations  may  now  be  drawn  regarding  the 
nature  of  the  normal  force  as  a  function  of  the  remaining  system  parameters 
and  variables: 

(a)  As  the  rotor  speed  increases: 

(i)  The  normal  force,  n,  first  becomes  positive  when: 

U62)2  =  [(1  -  82)2  +  (25S)2] 

This  simply  states  what  was  earlier  determined;  that  rub  initiates  when  the 

magnitude  of  the  centrifugal  force  due  to  mass  eccentricity  equals  the 

rotor's  impedance  multiplied  by  the  radial  clearance. 

2 

(ii)  The  value  of  ( 8  -  1)  increases;  which  tends  to  increase  the 

normal  force. 

2 

(iii)  The  value  of  the  centrifugal  force,  (ca  )  increases;  which  will 
increase  tne  normal  force  if  tne  phase  angle,  $,  is  less  than  90°;  and  will 
decrease  the  normal  force  if  if  exceeds  90°. 
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The  multiplicity  of  rotor  speeds  (roots)  which  make  the  radical  term  in 
Equation  (15)  identically  equal  to  zero  lead  to  three  possible  sequences  of 
rotor  modes  of  motion  occurring  within  the  rub  regime: 

(i)  If  all  the  roots  to  Equation  (15)  occur  at  speeds  below  the  speed 
at  which  rub  initiates,  then  a  real  solution  to  the  rub  equation  exists  for 
all  rotor  speeds.  This  does  not  preclude  the  possibility  that  rub  may  cease 
(i.e.,  n  <  0),  but  rather  indicates  that  the  rotor  dynamics  will  be  governed 
by  the  relationships  previously  developed  for  synchronous  whirl  (i.e., 
Equations  (4)  and  (12)). 

(ii)  If  one  of  the  roots  to  Equation  (15)  occurs  at  a  speed  greater  than 
the  rub  speed,  then  the  rub  equation  has  no  solution  for  rotor  speeds 
greater  than  this  limiting  speed  8j  .  It  is  possible  that  8  could  occur  at 
a  speed  for  which  rub  is  not  in  effect  (i.e.,  the  normal  force  is  less  than 
zero  according  to  Equation  (12)).  This  condition  shall  be  discussed  in 
greater  detail  in  a  later  section. 

(iii)  If  two  of  the  roots  to  Equation  (15)  occur  at  a  speed  greater  than 
the  rub  speed,  then  there  is  a  limited  regime  within  which  synchronous  rub 
is  not  in  equilibrium.  For  all  rotor  speeds  less  than  the  first  root  and 
greater  than  the  second  root,  the  rotor  dynamics  may  be  described  by 
Equations  (4)  or  (12).  Otherwise  synchronous  rub  is  not  a  possible  mode  of 
rotor  motion. 

Throughout  the  speed  regime  in  which  synchronous  rub  is  not  a  solution, 
three  modes  of  rotor  motion  are  possible  -  either  a  forward  synchronous 
whirl  without  rub,  a  partial  rub  or  a  full  annular  rub  reverse  whirl.  These 
cases  shall  be  discussed  in  greater  detail  in  the  following  section.  In 
order  to  provide  clarity  on  the  interrelationships  of  these  various  modes, 
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the  outline  in  Figure  10  provides  a  sequencing  of  the  possible  system 
dynamics  as  a  function  of  rotor  speed.  It  is  informative  to  demonstrate 
these  various  modes  by  comparing  the  amplitude  of  displacement  and  normal 
force  experienced  by  the  rotor  as  functions  of  rotor  speed  as  in  Figures  11 
through  111. 

(i)  No  rub 

Tne  reference  case  is  the  simple  condition  that  rub  never  occurs.  Thi 
is  true  for  systems  where  e  <  2x,.  As  can  be  seen  in  Figure  11,  the  radial 
clearance  is  greater  than  the  maximum  displacement  experienced  by  the  rotor 
and  so  the  normal  force  is  zero  for  all  speeds. 

(ii)  Limited  synchronous  rub 

As  depicted  in  Figure  12,  rotor  rub  initiates  at  some  speed  8rub  <  1.0. 
Further  increases  in  rotor  speed  drive  the  rotor  hard  against  the  stator 
while  the  normal  force  steadily  increases.  The  large  friction  and  damping 
forces  which  also  develop  ultimately  disrupt  the  synchronous  rub.  Further 
increases  in  rotor  speed  do  not  lead  to  a  restoration  of  rotor  rub  but 
rather  follow  the  course  of  action  for  the  free  rotor. 

(iii)  Synchronous  rub-limited  discontinuity 

Similar  to  the  limited  synchronous  rub  case,  rub  initiates  at  a  speed 
BrUb  <  and  ceases  at  some  speed,  8^  .  At  Bj  ,  rotor  precession  jumps 

down  to  the  free-rotor  synchronous  whirl.  As  shown  in  Figure  13,  at  a  yet 
higher  speed,  dynamic  conditions  once  again  become  favorable  for  a  stable 
synchronous  rub.  At  this  speed  though,  some  perturbance  is  required  to 


bring  the  rotor  into  contact  with  the  stator  and  thus  initiate  the  rub.  It 
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must  be  emphasized  that  reverse  whirl  may  also  be  stable  at  this  speed 
(dependent  upon  the  slip  condition),  thus  making  it  difficult  to  predict  the 
exact  motion  which  would  ensue. 

(iv)  Continuous  synchronous  rub 

As  in  the  limited  rub  case,  rotor  rub  initiates  at  some  rotor  speed, 

6  <  1.0.  As  rotor  speed  is  further  increased,  the  rub  reaction  force 

rub 

also  steadily  increases.  In  this  case  though,  the  damping  and  friction 
forces  are  not  significant  enough  to  disrupt  the  equilibrium  of  the 
synchronous  rub  (see  Figure  1^). 
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Analytically,  regimes  where  the  various  modes  dictate  the  behavior  of 
the  rotor  may  be  defined  and  mapped  onto  a  family  of  curves  using  variables 
thus  far  defined.  Figure  15  depicts  the  impact  on  rotor  motion  of  varying 
the  coefficient  of  friction,  unbalance  (or  inversely  the  radial  clearance), 
or  speed  ratio  for  a  fixed  value  of  rotor  damping.  For  this  particular  case 
of  t  =  .10,  region  (i)  of  the  plot  depicts  the  condition  for  which  no  rub 
occurs  (independent  of  the  rotor  speed)  and  is  described  by  the  inequality 
of  e  <  2 r.  For  a  particular  value  of  e  >  2£,  rub  first  initiates  at  the 
speed  defined  by  the  previously  derived  relationship: 

_  (2£2:  D  ♦  X252:  i)2+_(e2-  12 

rUb'  (e2-1) 

This  value  provides  the  lower  bound  to  region  (ii)  in  Figure  15.  Region 
(ii)  depicts  the  speed  regime  wherein  synchronous  rub  is  in  equilibrium.  If 
we  consider  the  specific  case  where  the  coefficient  of  friction,  y  =  0.25, 
then  region  (iii)  depicts  the  speed  regime  where  synchronous  rub  has  fallen 
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FIGURE  10.  ROTOR  MODES  OF  MOTION  OUTLINE 
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out  of  equilibrium.  The  lower  bound  to  this  region  is  defined  by  the 
inequality  provided  in  Equation  (15): 


(e£T) 


2,2  _  LyU_z_I?)_t_i25B)a?  < 

2 

O  +  U  ) 


Synchronous  rub  never  re-establishes  equilibrium  in  rotational  systems  with 
a  value  of  e  described  by  the  following  inequality: 


0  +  v  ) 


Otherwise,  the  speed  at  which  synchronous  rub  recurs  is  depicted  by  the 
curve  labelled  "y  =  0.25."  This  curve  is  generated  similarly  to  Equation 


(15)  above,  but  now  represents  the  case: 

ue2>2  -  :.(2S6212  s  0 

(1  +  U  ) 

Lastly,  for  any  system  described  by  ?  ■=  0.10  and  y  =  0.25;  if 
22  2  1/2 

e  2  (5  +  y  )/[y(1  +  y  )  ],  then  synchronous  rub  is  a  solution  to  the  rotor 

motion  for  all  speeds,  6  >  6  .  A  plot  such  as  in  Figure  15  thus  clearly 

outlines  rotor  motion  as  a  function  of  speed  for  any  system  described  by  the 
variables  y,  z,  and  e.  Specifically,  each  of  the  four  cases  described  in 
this  section  are  delimited  according  to  the  system's  value  of  c: 

(i)  No  rub:  e  <  2c; 

(ii)  Limited  synchronous  rub: 


2$  <  e  S 


(1+W2) 


(iii)  Synchronous  rub  -  limited  discontinuity: 


<  e  < 


( 1  +  y2 ) 


2  2 
y( 1 +y2) 


i).  DYNAMICS  OF  REVERSE  WHIRL 


Reverse  whirl  is  the  state  of  rotor  dynamics  wherein  the  direction  of 
precession  of  the  rotor  center  about  the  bearing  axis  opposes  the  direction 
of  rotation  of  the  rotor  about  its  own  center.  It  was  shown  that  in  the 
absence  of  rub,  an  unbalance  in  the  rotor's  mass  distribution  results  in  a 
forward  synchronous  whirl  of  the  rotor.  The  presence  of  damping  resists 
this  forward  precession  and  thus  creates  a  phase  lag  between  the  rotor 
center  as  it  precesses  and  the  mass  center  as  it  rotates.  This  damping 
force  though,  is  linearly  related  to  the  whirl  speed  and  so  can  not  by 
itself  destabilize  the  precessional  motion  of  the  rotor.  When  rub  occurs,  a 
frictional  force  at  the  interface  between  the  rotor  and  stator  develops 
which  combines  with  tne  damping  force  to  oppose  the  rotor's  forward 
precession.  Unlike  the  damping  force  though,  the  frictional  force  is  a  non¬ 
linear  function  of  rotor  speed.  If  the  additional  force  created  by  the 
frictional  force  opposing  the  rotor's  forward  whirl  can  not  be  balanced  by 
the  tangential  component  of  the  centrifugal  force  associated  with  the 
rotor's  mass  unbalance,  then  the  synchronous  whirl  assumption  breaks  down. 

The  dynamics  of  the  reverse  whirl  problem  are  far  more  complex  than  for 
forward  synchronous  whirl.  Since  the  whirl  speed,  fi,  is  now  independent  of 
the  rotational  speed,  w,  it  must  be  introduced  as  a  new  variable  in  the 
analysis.  This  then  creates  the  additional  complexity  of  having  a  phase 
angle  which  is  now  a  function  of  time  (more  specifically,  a  function  of 
(w  -  Jl) t) .  Since  the  phase  angle  varies  regularly  with  time  though  (i.e., 
constant  frequency),  a  simplifying  assumption  is  that  the  time  averaged,  or 
net  position  of  the  mass  center  of  the  rotor  is  located  at  the  centroid  of 
the  rotor  (i.e.,  6=0;.  Furthermore,  since  the  whirl  and  rotational 
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directions  are  opposed,  then  the  direction  of  relative  motion  at  the  point 
of  contact  (designated  point  A)  between  rotor  and  stator  must  be  determined. 
The  velocity  of  point  A  may  be  written: 


V  =  V  +  V. 
c  Ac 


where  |V  |  =  -  =  linear  velocity  of  the  rotor  center 


| V Ac  1  =  aw  =  linear  velocity  of  point  A  relative  to  the  rotor  center 
where  a  =  rotor  diameter 

By  assuming  (au)  >  (r  ft) ,  then  the  linear  velocity  of  the  rotor  at  A  may  be 
assumed  in  the  direction  of  the  rotor's  rotation  and  the  friction  force  is 
thus  aligned  with  the  direction  of  whirl. 

The  free  body  diagram  may  now  be  analyzed  to  determine  the  conditions 
whereby  reverse  whirl  is  possible.  Figure  16  displays  the  forces  at  work 
and  their  assumed  orientation  for  reverse  whirl.  By  considering  the  balance 
of  forces  in  the  tangential  and  radial  directions,  two  independent  equations 
are  developed.  In  the  tangential  direction,  the  friction  force  is  balanced 
solely  by  the  damping  force: 

uN  =  D  (r  ft) 
r  c 

By  defining  p  =  ft/oo  and  substituting  non-dimensional  variables: 

un  =  2  c  B  p  (16) 

In  the  radial  direction,  the  centrifugal  force  of  the  rotor  is  balanced  by 
the  elastic  force  of  the  rotor  plus  the  normal  force  between  the  rotor  and 
stator: 

Mrft2  =  N  +  Kr 
r  c  r  c 

In  non-dimensional  form  this  becomes: 

p2B2  =  n  +  1  (17) 

The  rotor  rotational  speed,  w,  and  whirl  speed,  ft,  may  be  kinematically 


*'/'//  V  .V.V  ~ *-•  V 

m V.*' * 
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related  through  the  slip  condition  that  exists  at  the  rotor-stator 

interface.  A  no-slip  condition  would  yield  the  kinematic  relationship 

common  to  epicyclic  gear  trains: 

-  (a/r  )w 
c 

If  there  is  slip  at  the  interface,  the  whirl  speed  is  reduced  accordingly. 
The  slip  ratio,  s,  is  thus  defined  as  the  ratio  of  the  (rotor  whirl  speed)- 
to-the-  (rotor  whirl  speed  without  slip);  i.e., 


(a/r  )oj 
c 


(17a) 


Equations  (16)  and  (17)  may  thus  be  rewritten: 

pn  =  2tB[(a/r  )s]  (16a) 

[s(a/rc)]262  =  n  +  1  (17a) 

By  substituting  (17a)  into  (16a),  the  critical  condition  for  reverse 
whirl  is  determined: 

u{[(a/rc)2s2]e2  -  1)  =  2;6[(a/rc)s] 

The  above  condition  may  be  resolved  into  a  quadratic  equation  in  B: 

y[(a/rc)s]262  -  2;[(a/rc)s]6  -  p  »  0  (18) 

Solving  for  6  yields  the  critical  speed  above  which  reverse  whirl  is 
possible: 

2  2  1/2 
r  .  C+U  ♦  U  ) 

-  [(");][ - ]  (19) 

rev  as  p 

The  lone  unknown  variable  in  the  above  equation  is  the  slip  ratio,  s, 
which  may  in  fact  be  difficult  to  determine.  The  limiting  case  arises  of 
course  when  no  slip  is  present  (i.e.,  s  =  1.0).  If  it  is  assumed  that  the 
damping  and  friction  are  of  the  same  order  of  magnitude,  then  8rey  is  of  the 
order  of  r  /a,  which  in  general  is  much  less  than  unity.  The  limiting  case 
thus  states  that  reverse  whirl  is  possible  for  essentially  all  rotor  speeds 


_  J 
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of  interest  (i.e.,  within  the  rub  regime).  It  is  not  until  a  significant 
degree  of  slip  is  present  (i.e.,  s  is  of  the  order  of  r^/a)  that  the  minimum 
possible  speed  for  reverse  whirl  may  become  restrictive  on  the  behavior  of  a 
rotor  undergoing  full  annular  rub. 

Since  the  rotor  first  assumes  a  forward  synchronous  precession  though, 
it  will  continue  with  this  motion  until  either  it  becomes  unstable  or  is 
perturbed  into  another  stable  mode  of  motion.  As  previously  developed,  the 
forward  synchronous  whirl  without  rub  is  continuously  stable  as  long  as  the 
displacement  of  the  free  rotor  is  less  than  the  constrained  rotor's  radial 
clearance.  The  only  mode  where  synchronous  whirl  may  be  destabilized  then, 
is  during  a  synchronous  rub.  This  has  been  analytically  developed  and 
occurs  if  the  reaction  force  between  the  rotor  and  stator  becomes  imaginary 
or  if  the  equilibrium  state  for  the  rotor-stator  is  inherently  unstable. 

In  this  situation,  the  equations  developed  based  on  synchronous  whirl  are  no 
longer  applicable. 

A  qualitative  analysis  of  the  transient  state  following  destabilization 
of  the  synchronous  rub  proves  quite  meaningful.  It  may  first  be  reasonably 
presumed  that  the  rotor  will  assume  either  a  partial  rub  or  a  forward 
subsynchronous  whirl  following  the  destabilization.  Considering  the  latter 
mode  of  motion,  the  impact  of  the  non-synchronicity  on  the  rotor's  dynamics 
is  primarily  that  the  phase  angle,  4> ,  is  now  a  time  dependent  variable.  The 
centrifugal  force  which  has  been  pressing  the  rotor  against  the  stator 
throughout  the  full  annular  rub  now  has  a  fluctuating  component.  As  the 
phase  angle  approaches  180°,  the  centrifugal  force  acting  on  the  rotor 
center  due  to  the  mass  unbalance  is  actually  joining  with  the  spring  force 
to  drive  the  rotor  center  in  towards  the  bearing  axis.  The  normal  force 
rapidly  drops  to  zero  and  contact  is  broken  between  the  rotor  and  stator. 


This  may  be  somewhat  anticipated  by  observing  the  near-vertical  negative 
slope  of  the  corresponding  normal  force  plot  at  the  critical  equilibrium 
speed  (see  Figure  17).  Now  that  rub  is  no  longer  present,  there  is  no  force 
to  drive  a  reverse  whirl,  and  so  the  rotor  will  jump  down  to  the  stable  mode 
of  forward  synchronous  whirl  without  rub.  Perturbing  the  rotor  back  into 
contact  at  this  speed  (or  higher)  may  now  provoke  a  reverse  whirl.  In  the 
previous  section,  the  possibility  for  a  restabilization  of  synchronous  rub 
was  demonstrated.  Under  these  conditions,  either  synchronous  rub  or  reverse 
whirl  may  result  when  the  rotor  is  perturbed  back  into  contact  with  the 
stator.  In  general,  no  smooth  transition  can  occur  from  a  synchronous  whirl 
to  a  reverse  whirl  since  this  would  mean  transiting  through  a  regime  of  zero 
whirl  speed.  The  whirl  speed  though  would  not  drop  below  that  necessary  to 
maintain  rub-contact;  for  when  contact  is  lost,  the  rotor  either  assumes  a 
partial  rub  whirl  or  resumes  a  synchronous  whirl  without  rub. 

The  behavior  of  the  rotor  when  transiting  from  a  higher  to  a  lower 
speed  differs  markedly  from  the  behavior  of  the  rotor  as  it  is  increasing 
speed.  In  the  absence  of  some  outside  perturbation,  when  reducing  speed 
from  a  situation  of  synchronous  whirl  without  rub  at  a  speed,  B  >  1.0;  no 
jump  phenomenon  in  the  rotor  displacement  curve  would  occur,  but  rather  rub 
would  occur  only  within  that  region  of  Figure  (5b)  where  the  free-rotor's 
displacement  exceeds  the  rotor  radial  clearance. 
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ROTOR  SPEED 


FIGURE  17.  NORMAL  FORCE  PLOT  O.Wt*d  Synchronous  Rub) 
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5.  ROTOR  RUB  WITH  A  NON-RIGID  STATOR 


The  assumption  of  an  ideally-rigid  stator  enclosing  the  rotating 
element  in  this  analysis  is  not  feasible,  but  depending  on  the  stator's 
dynamic  characteristics,  such  may  be  an  entirely  satisfactory  assumption  for 
the  analysis.  In  this  section,  the  stator  is  assigned  a  finite  mass,  radial 
stiffness  and  damping  coefficient.  For  ease  of  analysis,  these  characteris¬ 
tics  are  assumed  to  be  axisymmetric.  Additionally,  the  stator  is  assumed  to 
be  torsionally  rigid.  It  is  often  instrumental  to  investigate  the 
translational  system  analogous  to  the  rotating  system  under  analysis.  This 
is  first  accomplished  to  provide  insight  to  the  complex  equations  developed 
in  the  final  analysis. 

The  basic  precept  in  analysing  the  coupled  rotor-stator  system  is  that 
once  the  state  of  continuous  contact  is  interrupted,  the  synchronous  full 
annular  rub  condition  will  rapidly  deteriorate  to  a  partial-  or  no-rub 
state;  or  if  sufficiently  perturbed,  to  a  reverse  whirl.  The  translational 
system  designed  to  model  the  rotor-stator  is  presented  in  Figures  1 8 ( a )  and 
(b).  The  rotor  and  stator  are  assigned  their  respective  masses,  a  linear 
spring  and  damper  to  characterize  their  dynamic  characteristics,  and  are 
restricted  to  motion  along  the  horizontal  axis.  In  the  rotational  system, 
steady  motion  of  the  combined  rotor-stator  may  be  depicted  by  an  arc  which 
traces  the  path  of  the  rotor-stator  center  of  mass.  The  actual  force 
experienced  by  the  combined  mass  would  have  a  frequency  equal  to  the  rotor's 
precessional  frequency,  and  a  magnitude  which  depends  upon  the  centrifugal 
force  acting  on  the  rotor  and  the  dynamic  characteristics  of  both  the  rotor 
and  stator.  In  the  translational  system  a  static  pre-load,  F  ,  is  presumed 
to  displace  the  rotor  mass  across  the  radial  clearance,  r  ,  and  then  further 


E 


displaces  the  combined  rotor-stator  mass  to  a  static  equilibrium  position 
about  which  the  system  vibrates.  A  sinusoidal  force  then  drives  the  system 
at  the  'synchronous'  speed,  w  (analogous  to  the  rotor's  rotational  speed), 
with  a  magnitude  F  . 

The  free  body  diagrams  for  the  rotor  and  stator  under  the  static  pre¬ 
load  are  displayed  in  Figure  18(c). 

The  force  which  acts  at  the  interface  is  given  by  P  .  Essentially, 
only  the  stiffness  of  the  rotor  and  stator  resist  the  static  load.  The 
static  equation  for  the  rotor  and  stator  are  given  by  Equations  (19)  and 
(20)  respectively: 


F  -  P  =  k, x. 
o  o  11 


Po  =  k2X2 


Additionally,  x1  and  x2  are  related  by  the  kinematic  condition: 


1 


r  +  x., 
c  2 


Therefore,  a  single  equation  may  be  written: 

F  -  P  -  k.(r  +  P  /k.) 
o  o  1c  o  2 

Solving  for  the  static  load  at  the  interface  yields: 

P  =  (F  -  k  r  )k  /(k  +  k  ) 

o  o  1  c  2  1  2 


(19) 

(20) 


(21) 


The  equation  of  motion  for  the  system’s  dynamic  condition  is  similarly 
determined  by  analyzing  the  forces  acting  on  the  rotor  and  stator.  The 
dynamic  condition  free  body  diagrams  are  provided  in  Figure  18(d). 


Applying  Newton's  2nd  Law,  the  balance  of  forces  may  be  written: 
rotor: 


-k  x  -  D  x  -  P  -  P  e^wt  +  F  +  F  e^ut  =  m  x 
1  1  1  1  o  1  o  1  11 


(22) 


stator: 


— k  x  —  D  x  +  P  +  P  p 
2  2  2  2  o  1 


jut 


m2x2 


where  P1eJa)t  represents  the  dynamic  reaction  force  at  the  rotor-stator 


(23) 


STATOR 


interface. 


By  assuming  the  coupled  system's  response  is  a  harmonic  motion  of  magnitude 


A,  and  recalling  the  kinematic  relationship  between  x1  and  x^,  the  rotor's 


and  stator's  displacement  variables  may  be  rewritten: 


Jut. 


x.  =  (r  +  P  /k  +  Ae  ) 
1  c  o  2 


x.  =  (P  /k  +  AeJa)t) 
2  o  2 


Making  these  substitutions  into  Equations  (22)  and  (23)  results  in  the 

following  equations  of  motion: 

rotor: 


-k1  (  rc  +  PQ/f<2  +  Ae';’wt)  -  D1  ( ju)  Ae''t,'r-Po  -  P1  e 


jwt+ 


F  +  F  e 
o  1 


J  Ut  rr,  ,  ^  JUt 
—  “175  ^  &0J  <2 


stator: 


-k. (P  /k  +  Ae'-'“t)  -  D_( jw) Ae''“t+  P  +  P  e^Uit  =  m.,Au2e^U)t 


2'  o  2  2W  '  o  1  2 

These  relationships  are  valid  only  if  the  rotor  and  stator  elements  maintain 


continuous  contact.  Eliminating  the  static  equilibrium  relationships  from 
the  above  equations  and  solving  each  equation  for  the  dynamic  reaction  at 


the  interface  yields  equivalent  relationships: 


P1  =  F1  +  A ( — k  ^  -  juD  +  m  u  ) 


A(k2  +  juD2  -  m2u  ) 


These  two  relationships  may  be  reduced  to  a  single  equation  and  solved 
for  the  magnitude  of  the  vibrational  response: 


A  = - r— 1 -  (24) 

[  ( k  1  +  k2)  -  u  '-J+  m2>]  +  ju(D1+  D  ) 

The  criterion  that  'synchronous  rub'  (as  portrayed  by  this  model) 
remains  in  equilibrium  is  evidenced  through  the  model  by  a  positive 
reaction  force  at  the  interface,  i.e.,  that  the  magnitude  of  the  dynamic 


reaction  force  be  less  than  the  static  reaction  force: 


-  «... 


disregarded  since  only  the  magnitudes  involved  affect  the  "equilibrium 
criterion".  Rearranging  terms  results  in  the  final  criterion  dictating  the 
equilibrium  of  synchronous  rub: 


(F  -  kiO 
o  1  c 


[(1  “  ~%)2  +  (2c  --)2] 

2  0  “h 

>  N 


2  .  ,1/2 

[(1  "  ~—2)2  +  (2^2  -”)2] 


(26a) 


The  inter-relationships  of  the  various  terms:  | F ^ )  relative  to 
(F  -  ^r  ),  and  {[1  -  (w2/^2  )]2  +  [25q( w/u^ ) ] 2} 1 /2  relative  to 
222  21/2 

{[1  -  U  /'o)N2  )]  +  [2c, ( w/ui^)]  }  outline  various  regimes  in  the 

plot  of  (force  vs.  frequency)  where  synchronous  rub  is  in  equilibrium.  Four 
fundamental  cases  are  depicted  graphically  in  Figures  (19)  through  (22): 
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Throughout  tnese  figures,  the  shaded  region  indicates  the  speed  regime 

wherein  synchronous  rub  does  not  satisfy  dynamic  equilibrium  in  accordance 

with  the  condition  provided  by  Equation  (26a).  The  artificialities  of  the 

translational  model  are  evidenced  by  the  awkward  low  frequency  behavior  of 

the  graphs.  In  particular,  there  is  no  clear  correlation  between  the  terms 

(F  -  k.r  )  and  F, ,  and  similar  functions  in  the  rotational  system.  A 
o  1  c  1 

single  most  prominent  observation  in  all  four  cases  though,  is  that 
synchronous  rub  is  not  in  equilibrium  in  the  vicinity  of  the  combined  rotor- 
stator  natural  frequency.  This  supports  earlier  statements  that  it  is  the 
characteristics  of  the  combined  rotor-stator  which  govern  the  stability  of 
synchronous  rub. 


Combining  Equations  (22),  (23)  and  (24),  this  criterion  becomes: 


|F.,-  A[(k1~  id  m1  )  +  JojD1  ]  |  <  (Fq-  k1rc)k2/(k1+  k2) 


(25) 


|  A[ ( k2~  u)  m2)  +  J(dD2]  |  <  (Fq-  k1rQ)k2/(k1+  k2) 


(26) 


Substituting  for  the  variable  A  using  Equation  (24),  this  condition  becomes 

F.  [ ( k  -  oj2m  )  +  ju)D  ]  k 

|F - 1 - 1  __1 - 1 -  |  <  (F  k  )(  i 

t(k1+  k2)  -  a  (m1+  m2)]  +  Dg)  o  I  c  k,  k2 

It  is  convenient  to  employ  the  following  definitions  when  referring  to  the 
rotor,  stator  and  combined  rotor-stator: 

1  /  2 

=  [(k  +  k2 ) / (m  +  m2 ) ] 


“N1  =  (k1/mi) 


1/2 


“N2  =  (k2/m2) 


1/2 


40  “  ^D1+  D2)/2[(k1+  k2)(m1+  m2)] 


1/2 


C1  =  D1/2[k1m1  ] 


1/2 


C2  “  D2/2[k2m2] 


1/2 


The  condition  that  "full  annular  synchronous  rub"  be  a  solution  to  system 
motion  now  becomes: 


[(1  -  “-)  ♦  j(2c2  )] 

If  k2  _  _  _  ^2 _  WN2_  ,  k2 

lFl'  k  +  k  2  I  (Fo  kirc)  k'+'k" 

[(1  --U-2)  ♦  i(2,0^-)] 

“N  N 


This  inequality  may  be  further  simplified  by  eliminating  common  factors  and 
by  writing  the  two  complex  dimensionless  impedances  in  terms  of  their 
magnitude  and  phase.  The  phase  relationships  which  evolve  may  be 
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The  more  exact  analysis  of  the  rotor-stator  dynamics  requires  the 
assignment  of  finite  mass,  radial  stiffness  and  damping  characteristics  to 
the  stator  in  Figure  2.  As  earlier  stated,  these  characteristics  are 
assumed  axisymmetric .  The  reference  position  for  both  the  rotor  and  stator 
displacements  remains  the  imaginary  bearing  axis,  B,  which  marks  the  rest 
position  for  the  rotor  and  stator.  Figure  23  depicts  the  rotor  in  contact 
with  the  stator;  with  the  free  body  diagrams  for  each  element  provided  in 
Figures  24(a)  and  (b).  The  equations  of  motion  for  the  rotor  and  stator  are 
developed  according  to  the  balance  of  forces  in  the  free  body  diagrams.  The 
forces  acting  on  the  rotor  are  identical  in  nature  to  the  forces  present  in 
the  rigid  stator  case.  The  balance  of  forces  for  the  rotor  is  written: 


-K  2  -D  2  -R(t)  =  M  (2  ) 

r  c  r  c  r  g 


(27) 


The  reaction  force  which  acts  at  the  rotor-stator  interface  though,  is  now  a 
variable  which  depends  on  both  the  rotor's  and  stator's  dynamic 
characteristics.  Since  the  stator  is  non-rigid,  the  reaction  force  drives 
the  stator  out  of  line  with  the  bearing  center-rotor  center  axis.  The 
resultant  phase  relationship  between  the  rotor  displacement  and  the  rub 
location  may  be  described  by  the  angle  X,  which  is  a  function  of  the  rotor 
displacement,  stator  displacement  and  reaction  force  lines  of  action. 

Tne  reaction  force  may  thus  be  broken  down  into  its  normal  and  tangential 
(friction)  components  with  a  periodic  time  dependence  characterized  by  the 
synchronous  speed,  w,  and  the  phase  angle  X: 


R  ( t)  =  (N  -  j  yN )  e 


j  ( wt  +  X ) 


or 


where 


R ( t)  =  N ( 1  +  y  )  e 


ip  =  tan  ^  y 


+  2.1/2  jUt  +  X  +  *) 
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as 


Additionally, 

expressed: 

z 

g 


the  displacement  of  the  rotor's  center  of  gravity  may  be 


5ejUt+6) 


6  is  the  phase  angle  which  accounts  for  the  lag  between  the  rotor  center 
and  its  forcing  function.  This  forcing  function  now  comprises  both  the 
centrifugal  force  due  to  unbalance,  F  ,  and  the  reaction  force,  R(t).  The 
resulting  expression  for  6  is  thus  a  function  of  three  phase  angles:  <J>, 
the  phase  angle  due  to  rotor  damping;  <jj,  the  phase  angle  describing  the 
line  of  action  of  the  reaction  force  at  the  rub  location;  and  \,  the  phase 
difference  between  the  rub  location  and  the  rotor  center's  displacement: 

6  =  tan-1[u)Dr+  N  ( 1  +  y2) 1  /2Sin(  X ) ]/ [  ( K  -  u>2M  )  +  N  ( 1  +  y2) 1  /2cos(  X )  ] 

Making  these  substitutions,  the  force  balance  equation  for  the  rotor  may 
now  be  written: 


-K  z  -  D  z  -  N(1  +  y2)1,/2eJ(tjt  +  A  +  ^)  »  M  -~(z  +  6ejUt+0))  (28) 

r  c  r  c  r  2  c 

at 

Equation  (28)  may  be  solved  for  the  displacement  of  the  rotor  center: 
uj2M  6ejUt+6)-  N(1  +  u2)1/2ejwt+X  +  4,) 


L(Kr-  u>  Mp)  ♦  (juDp)] 


Breaking  the  rotor's  impedance  down  to  its  magnitude  and  phase  yields: 

2m  .  jUt+6)  NM.  2.1/  2  j  (  ut+  A+iJ>) 
oj  M  6e  -  N  ( 1  +  y  )  e 


z  = 
c 


2  2  21/2 
C(Kr-  a\r  +  (ud r)  r 


(29) 


Dividing  through  by  the  term  rce^ut  and  substituting  non-dimensional 
variables,  the  rotor  displacement  may  be  finally  written: 


-► 


z 

c 


r 

c 


U)t 


[(1-S2)2  +  (2tB)2]1/2 


(29a) 


Whereas  the  centrifugal  force  due  to  unbalance  drives  the  rotor’s 


motion,  it  is  the  synchronous  rub  force  (i.e.,  the  reaction  force,  R(t)) 
which  acts  as  the  input  driving  force  for  the  stator's  motion.  The  forces 
which  respond  to  this  driving  force  are  the  stator's  inertial,  elastic  and 
damping  responses.  The  force  balance  for  the  stator  is  derived  from  the 
free  body  diagram: 


-K  z  -Dz-R(t)=Mz 
s  s  s  s  s  s 


By  assuming  a  synchronous  whirl  response  by  the  stator,  z  =  |z  |e 

s  s 


j( wt+n) 


and  solving  for  zs,  the  stator's  equation  of  motion  is  determined: 


Now  dividing  through  by  the  term  r  e  and  factoring  out  the  ratio  (K  /K  ), 

c  r  s 


the  non-dimensional  relation  for  the  stator's  displacement  becomes: 


r  e 
c 


s _ 

jut 


,  r.  2,1/2  jU+X-t) 
(--)n(1+p  )  eJ 

s 


2  2  2  1/2 
[(1-BS  )  ♦  (2?s6s)‘f]1/^ 


(31  ) 


Two  variables  of  interest  in  this  analysis  are  the  rotor's  displacement 
and  tne  reaction  force  at  the  rotor-stator  interface.  To  determine  the 


rotor's  displacement  z ,  Equations  (27)  and  (30)  are  utilized.  The  reaction 


force  acting  upon  the  stator  is  equal  in  magnitude  but  opposite  in  direction 
to  that  acting  upon  the  rotor,  so  Equations  (27)  and  (30)  may  be  combined 
with  the  force  term  R(t)  eliminated: 


-K  z  -K  z  -D  z  -D  z  =  M  z  +  M  z 
rc  ssrcss  rg  ss 


(32) 


Figure  23(b)  provides  the  vector  triangle  relating  the  stator  displacement, 
rotor  displacement  and  radial  clearance.  The  kinematic  condition  which 
arises  from  the  geometry  presented  here  is  given  by: 

j(  tot+X) 


z  -  z  =  re" 
c  s  c 


The  implication  here  is  twofold: 

(i)  The  magnitude  of  the  rotor  displacement  relative  to  the  stator 
displacement  equals  the  radial  clearance  between  the  two  elements. 

(ii)  The  point  of  contact  between  the  two  elements  must  be  in  line  with 
the  centers  of  the  displaced  rotor  and  stator. 

Applying  this  kinematic  relationship  and  the  previously  derived  condition, 


♦  6eJUt+*) 


z  *>  z  +  5e  ,  results  in  the  following  equation  of  motion  for  the 

8  c 


rotor: 


(M  +  M  )z  +  (D  +  D  )z  +  (K  +  K  )z  =  M  ~~~(  6e^(  ut  +  <p) )  +  K  re 
rsc  rsc  rsc  r  Jk2  s  c 

dt 


j(wt+X) 


The  synchronous  whirl  solution  is  thus  given  by: 


[o>2M  6ej4l+  K  r  ejA]ejwt 
-►  r  S  C 

zc  2 

t(Kr+  Ks)  -  u  (Mr+  Mfi)]  ♦  j[to(Dr+  Dg)] 


The  time-independent,  non-dimensional  solution  for  the  rotor  displacement 
is  obtained  by  dividing  each  side  by  the  quantity  rceJut  and  factoring  out 
on  the  right  hand  side.  The  resulting  expression  is: 


J* 


e~jK  (33a) 
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{[(1-6")  -  --(1-6  ----)]  ♦  [2tB  +  2t6““]  } 

Kr  ^Ns  ^Ns 


where  <  =  tan 


ui(D  +  D  ) 
r  s 


(K  +  K  )  -  oo  (M  +  M  ) 

PS  PS 


This  equation  is  further  reduced  to  a  single  magnitude  and  phase  in  the 
following: 

r,  „2  Ks  .2  ,  n2  .  ±  Ks  .  ,  . 2-,  1  /2 

[  ( e8  003$  +  r-cosA)  +  (e6  sin$  +  --smA)  J 

Z  ft  ft  .  /  V 

-S----  - . . C - 5 . . -  (33b) 

rce  2  Ks  2  “Nr  2  V  2  1/2 

°  { [  ( 1  ~8  ;  +  -~(1-62  +  [2CB  +  2^6  -“]  ) 

CJN  s  “n  s 
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(e6  sin$  +  --sinA) 

-1  r 

where  v  =  tan  - r. - 

2  Ks 

(e6  cos$  +  --cosA) 
ft 

P 

The  reaction  force  is  more  difficult  to  solve  than  the  rotor  displace¬ 
ment,  but  provides  much  meaningful  information  regarding  the  rub  motion.  A 
positive  reaction  force  indicates  full  annular  synchronous  rub  is  in 
equilibrium.  Conversely,  a  negative  reaction  force  indicates  a  no-rub 
action,  while  an  imaginary  reaction  force  indicates  full  annular  synchronous 
rub  is  not  in  equilibrium  (i.e.,  not  a  solution  to  the  rub  equation).  In 


,ii 


order  to  solve  for  the  reaction  force  as  a  function  of  the  rotor  speed,  the 


kinematic  condition  relating  the  rotor  and  stator  motion  must  once  again  be 
applied: 


z  -  z 

_C _ s  _  jX 

jut 

r  e 
c 

The  rotor's  displacement  relative  to  the  stator  is  derived  from  Equations 
(29a)  and  (31 ) : 

K 

z  _  -1  ^  2  j(e-<t>)  ) 6  JT 

VtS  -  “T  2  -  n(1  +  ,2)1/2e^+^4  +  -ST7r-J 

r  eJut  R  1/2  R  U‘  R  1/2 

c  r  r  s 

where  R  =  [(1~62)2  +  (2?e)2] 
r 

and  Rs  =  [(1~8S2)2  +  (2cs8s)2] 

Applying  the  above  kinematic  condition,  this  equation  may  be  written: 

(-C)eJ(*-T> 

.  ,,  2  1/2[A^  +  _ls_„  jA  =  jX 

D  1/2  v  M  '  D  1/2  _  1/2  J 
K  K  n 

r  r  3 

Solving  for  the  phasor  yields: 

eJ1  -  (Rs,/£c62eJtS-*)j/ 

{(RrR3)1/2+  n(1  +  p2)1/2CRs1/2cos(4'-<())  +  (")Rr,1 /2cos(  ^-x)j  + 

S 

jn(  1  +  m2)1 /2[R  1 /2sin(i|)-4i)  +  (--)R  1 /2sin(  iJj-t)  ] } 

S  f\  r 
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By  defining 


n(l+u  ) 


sin(<{i-<j>)  +  ( — — )  R 


(R  R  )  +  n(1  +  u  )  “LR  cos(4)-$;  +  (n~)R  cos(<|r-x)J 

r  s  s  K  r 

s 


then  the  relationship  becomes: 


ejX  =  lRs1/2EB2ej(0"^o)}/ 


({(RrRs)1/2+  n(1  +  u2)1/2[Rs1/cos(4;-4>)  +  ( --)Rp1 /2os(  )  ] } 2 


,,  2,1/2r_  1/2  .  ,  .  ,  r,  D  1/2  .  ,,,2,1/2 

+ 1  n(  1  +  p  )  [R  sin(ip— 4>)  +  (j^_  )Rp  sin(ijj— t)]}  )  (3^ 

s 

Therefore,  the  phase  angle  may  be  defined  in  terms  of  known  variables: 

X  =  (6  -  <j>  -  o)  , 

witn  the  magnitude  of  Equation  (32)  being  equal  to  unity.  The  square  of 
the  magnitude  may  also  be  equated  to  unity,  yielding: 


R  (esV/(KRrR  )U2+  n(1  +  p2)1/2[R  ,/2cos(*-<»  +  (--)R  1 /2cosU-t  )  ] } 2 


K  '  r 
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+  {n(1  +  p2)1/2[R  1 /2sin(t|>-<i>)  +  (--)R  1/2sin(t~i)]}2)  =  1  (3^ 

S  l\  i 

s 

Equation  (3*ia)  is  resolved  into  a  quadratic  equation  in  the  variable  n  in 
the  following  expression: 

K  R  1 72  KR2 

n2(  1  +  p2) 1  /  2[  1  +  2( cos (  t— )  +  (--)  --]  + 

K  R  K  R 

s  s  s  s 

,  K  R  (R  -  e2Bi)) 

i  /p  r  r  v  ' 

+  2n[Rr  cos(  4^-0)  +  ( — )—  1/2cos(<|;-t)]  + - 2~W2~  =  ° 

K  R  ( 1  +  p  ) 

s  s 


Solving  this  equation  for  the  normal  force,  n,  yields  the  lengthy 


expression: 
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An  alternative  expression  is  derived  by  applying  the  following 
relationships: 


cos(  -p  -  <p)  = 


LLl:B_)_t_25gu] 
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t(,-8s  )  * 


(,.„2),/2  R  1/2 
s 


cos(t  -  0)  = 


(1-62)(1-6q2)  +  (2^6) (2t  3  ) 

S  SS 


R  1/2R  1/2 
r  s 


Making  these  substitutions: 


K  R 
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K  R 

±  {[(1-62+  2cBu)  -  ~  R-n-6s2+  25s6su)]2 
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(35) 
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It  should  be  noted  that  the  expression  above  reduces  to  the  earlier 

derived  result  for  rub  against  a  rigid  stator  (Equation  11))  by  employing 

the  corresponding  conditions: 

K  =  therefore  K  /K  =0 
s  r  s 


R  =1.0 
s 


The  impact  of  releasing  restrictions  on  the  stator's  dynamic  behavior 

is  particularly  complex  due  to  the  numerous  variables  involved.  By 

comparing  Equations  (11)  and  (35)  though,  qualitative  conclusions  may  be 

drawn.  There  are  three  dynamic  'elements'  significant  to  this  system:  the 

rotor,  the  stator  and  the  combined  rotor-stator .  The  character  of  the 

response  of  each  element  is  dependent  on  the  speed  of  excitation  and  the 

element's  natural  frequency.  By  considering  the  magnitude  of  a  generic 

22  21/2 

element's  mechanical  impedance,  [(K  -  w  M)  +  (wD)  ]  ,  then  it  is  apparent 

that  at  low  excitation  frequencies,  the  response  is  as  a  stiff  spring, 

characterized  by  the  value  of  K.  At  high  frequencies,  the  response  is  as  a 

massive  element;  characterized  by  the  value  of  M.  The  response  is  a  maximum 

in  the  vicinity  of  the  natural  frequency,  /K/M,  (specifically  at 
2  2  1/2 

w  =  [u^  -D  /2M]  ).  These  simple  principles  may  be  applied  to  help  explain 

the  dynamics  of  the  rotor-stator  system. 

Given  the  rotor's  and  stator's  characteristic  mass  and  stiffness; 
the  composite  element's  natural  frequen../,  u^,  occurs  between  and  : 


Three  simple  relationships  thus  exist  between  the  rotor  and  stator: 

(i)  the  stator's  natural  frequency  is  greater  than  the  rotor's  natural 
frequency 

(ii)  the  stator's  natural  frequency  is  less  than  the  rotor's  natural 
frequency 

(iii)  the  stator's  natural  frequency  is  of  the  order  of  the  rotor's 
natural  frequency 

Each  of  these  relationships  may  be  satisfied  by  several  conditions  relating 

the  elements'  mass  and  stiffness.  For  example,  is  unconditionally 

true  when  K  >  K  and  M  =  M  or  when  K  =  K  and  M  <  M  .  It  is  condi- 
s  r  -  s  r  s  r  -  s  r 

tionally  true  when  K  >  K  and  M  >  M  or  when  K  <  K  and  M  <  M  .  It  is 

s  r  -  s  r  s  r  -  s  r 

evident  that  a  large  number  of  cases  may  be  analyzed  based  on  the  relative 
masses  and  stiffnesses  alone.  This  is  exacerbated  when  variations  in 
unbalance,  friction  and  damping  are  investigated. 

If  wNs  is  sufficiently  greater  than  eo^ ,  then  when  rub  first  initiates, 
the  stator  response  resembles  the  infinitely  stiff  case;  with  the  resulting 
reaction  force  depending  directly  on  Kg.  In  general,  since  it  is  the 
stator's  stiffness  which  keeps  it  pressed  against  the  precessing  rotor  at 
speeds  below  ,  then  continuous  synchronous  rub  becomes  a  less  stable  mode 
of  motion  as  decreases.  As  rotor  speed  approaches  ,  the  rotor-stator 
system  becomes  very  compliant  and  continuous  synchronous  rub  will  cease 
(i.e.,  the  radical  term  in  Equation  (36)  becomes  negative).  As  the  rotor 
speed  further  increases  and  reaches  u  ,  the  stator  alone  is  extremely 


compliant.  Since  a  small  push  by  the  rotor  thus  brings  a  disproportionately 

large  response  by  the  stator,  then  continuous  contact  may  rapidly  diminish. 

(This  is  evident  in  Equation  (36)  since  the  value  of  Rs  tends  toward  zero  as 

Ss  approaches  1.0  and  the  negative  term  within  the  radical, 

2 

-(K^/K  )  (Rr/Rg),  becomes  dominant).  Throughout  this  speed  regime,  the 
rotor's  motion  (if  unperturbed)  is  identical  to  that  of  the  free  rotor.  At 
yet  higher  speeds,  the  stator's  mass  dominates  the  character  of  its 
response.  Using  limit  analysis,  the  stator' ^  behavior  once  again  approaches 
the  rigid  response  at  very  high  speeds  (6  >>  B  ).  The  criterion  for 
stability  of  synchronous  rub  at  very  high  speeds  is  identical  to  that  for  an 
infinitely  stiff  stator;  i.e.. 


if  e  > 


<l*u2)1/2 


then  synchronous  rub  is  stable  for  8  >>  B 


At  this  speed  though,  the  rotor  must  be  perturbed  to  initiate  the  rub. 


If  ojNs  and  are  equal,  then  the  stator's  response  resembles  the 


rotor's  response  for  any  particular  excitation  frequency.  Therefore,  the 
speed  at  which  the  combined  rotor-stator  is  most  compliant  coincides  with 
the  speed  at  which  the  free-rotor's  runout  is  greatest.  This  similarity  in 


response  tends  to  nullify  the  stator's  non-rigid  character  if  (K  /K  )  and 


(NWMg)  are  greater  than  1.0.  As  Ks  and  Mg  are  increased,  the  response 


rapidly  approaches  the  rigid  stator  case.  When  the  stator's  dynamic 
characteristics  equal  the  rotor's,  then  the  speed  regime  in  which  continuous 
synchronous  rub  is  a  solution  to  rotor  motion  is  identical  to  the  rigid 
stator  case  and  the  reaction  force  is  exactly  half  that  which  would  result 
if  the  stator  were  rigid.  As  the  stator  becomes  less  stiff  and  less 


massive,  synchronous  rub  becomes  less  stable.  In  the  limit,  for  low  values 
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of  (IWK^),  the  system  is  unable  to  sustain  a  continuous  rub  with  the 
precessing  rotor. 

If  the  stator's  natural  frequency  is  measurably  less  than  the  rotor's, 

then  the  speed  regime  in  which  the  stator  is  most  compliant  occurs  prior  to 

the  initiation  of  rub.  When  the  rotor  first  gains  contact,  the  stator's 

impedance  and  that  of  the  combined  rotor-stator  are  increasing  with  the 

square  of  the  speed.  The  response  of  the  system  rapidly  approaches  that  of 

the  rigid  stator  case  as  M  increases  and  thus  either  w,,  decreases  or,  for 
°  s  Ns 

a  constant  w,,  ,  K  increases.  If  in  fact  a  low  value  of  u>.,  is  due  to  low 
Ns’  s  Ns 

values  of  both  Ks  and  M  ,  then  may  be  approximated  by  and  the  rub 

response  is  as  described  previously  for  . 

The  impact  of  the  stator's  damping  characteristics  are  straightforward. 

Independent  of  the  natural  frequencies,  stiffnesses  or  masses  of  the  rotor 

and  stator;  any  increase  in  the  stator's  damping  coefficient  tends  to 

diminish  the  equilibrium  of  synchronous  rub.  Assuming  light  damping  for 

both  rotor  and  stator,  variations  in  stator  damping  has  little  effect  on  the 

magnitude  of  the  reaction  force  during  rub. 

In  summary,  if  the  stator's  natural  frequency  is  greater  than  the 

rotor's  natural  frequency,  then  synchronous  rub  is  not  in  equilibrium 

throughout  a  speed  regime  which  encompasses  the  stator's  and  the  combined 

rotor-stator's  natural  frequencies  (w^s  and  respectively).  This  speed 

regime  widens  for  a  given  value  of  as  the  stator  stiffness,  K  ,  and 

Ns  s 

stator  mass,  M  ,  are  decreased.  Similarly,  the  reaction  force  experienced 
during  continuous  rub  for  a  given  value  of  u^s  decreases  as  Kg  and  are 
decreased.  As  the  stator's  natural  frequency  is  reduced  towards  the  rotor's 
natural  frequency,  this  particular  speed  regime  where  synchronous  rub  is  not 
viable  diminishes  accordingly.  When  coNs  =  the  stability  of  synchronous 


rub  resembles  that  of  the  rotor  against  a  rigid  stator;  being  more  stable  if 

(K^/Kg)  (and  thus  (fWM  ))  is  less  than  unity  and  vice  versa.  As  may  be 

expected,  the  reaction  force  in  this  case  increases  as  K  and  M  increase. 

s  s 

As  the  values  of  to^  and  are  reduced  below  the  synchronous  rub  system 
response  becomes  as  stable  as  if  the  rotor  were  rubbing  against  a  rigid 
stator.  This  is  true  since  the  stator's  impedance  rapidly  increases  with 
the  square  of  the  excitation  frequency.  Regarding  stator  damping,  any 
increase  in  the  damping  coef f icienct,  D  ,  reduces  the  speed  regime  in  which 
synchronous  rub  may  occur.  In  general,  the  stator's  finite  mass,  stiffness 
and  damping  characteristics  serve  to  abate  the  negative  characteristics  of 
rotor  rub.  In  no  case  do  these  characteristics  intensify  the  reaction 
forces  between  the  two  elements  (when  compared  with  rub  against  a  rigid 
stator),  but  rather  they  tend  to  reduce  the  reaction  forces,  often 
sufficient  enough  to  create  a  dynamic  relationship  less  favorable  for  full 
annular  synchronous  rub  (and  too,  reverse  whirl). 

Appendix  A  presents  a  numerical  investigation  of  the  impact  of  varying 
the  parameters  discussed  in  this  section.  Through  the  efforts  of  Professor 
Udo  Fischer  of  the  University  of  Magdeburg,  Magdeburg,  East  Germany,  a 
Fortran  program  has  been  developed  which  clearly  demonstrates  the  time 
dependent  response  of  the  rotor  and  stator  in  accordance  with  the  equations 
of  motion  developed  here.  The  stability  of  the  synchronous  rub  precession 
in  those  speed  regimes  where  this  mode  of  rotor  motion  is  in  equilibrium  is 


clearly  demonstrated. 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

Mechanical  vibrations  in  rotating  machinery  result  from  any  mass 
unbalance  present  in  the  rotating  element.  Should  rotor  rub  occur  as  a 
result  of  these  vibrations,  the  integrity  of  the  system  may  be  significantly 
impaired.  Two  direct  solutions  to  this  problem  are  immediately  apparent: 

(i)  Reduce  the  mass  unbalance  which  drives  the  vibrations  in  the 


2 _ 


system,  or 

(ii)  Increase  the  radial  clearance  between  the  rotating  element  and  the 
stationary  element. 

Both  of  these  solutions  are  costly,  either  through  the  precision  required  in 
the  manufacturing  phase  or  efficiency  lost  throughout  the  operational  life 
of  the  system. 

The  dynamics  of  rotor  rub  are  quite  complex,  but  in  general  a  system 
may  be  resolved  to  one  of  four  behavioral  types  dependent  on  the  system's 
parameters: 

(i)  Synchronous  precession  without  rub 
(ii)  Partial  rotor  rub 

(iii)  Full  annular  rotor  rub 

(iv)  Reverse  whirl 

Most  rotating  machinery  are  designed  for  the  no-rub  condition.  Analyti¬ 
cally,  rub  may  be  anticipated  if  the  ratio  of  the  mass  eccentricity-to-the- 
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radial  clearance  exceeds  twice  the  damping  ratio  for  the  system.  Should  rub 
occur,  the  equilibrium  of  full  annular  rub  is  then  dependent  on  both  the 
rotor's  and  stator's  dynamic  characteristics,  the  system's  unbalance,  the 
radial  clearance  and  friction  experienced  between  the  two  contacting 
surfaces,  and  the  rotor's  speed. 

Initial  analysis  was  conducted  with  the  assumption  that  the  system's 
stator  is  infinitely  rigid.  Dependent  upon  the  system's  damping  ratio, 
displacement-to-clearance  ratio  and  frictional  coefficient;  the  rotor's 
response  as  a  function  of  speed  can  then  be  characterized  according  to  one 
of  the  following: 

(i)  Limited  synchronous  rub:  Synchronous  rub  is  in  equilibrium  up  to  a 
"critical  speed"  (unique  to  each  system)  above  which  equilibrium  is  never 
regained. 

(ii)  Synchronous  rub  -  limited  discontinuity:  Similar  to  the  limited 
rub  case,  equilibrium  is  lost  at  some  critical  speed,  but  is  then  regained 
at  a  higher  speed. 

(iii)  Continuous  synchronous  rub:  As  the  description  implies,  this  type 
system  is  able  to  find  an  equilibrium  position  for  all  rotor  speeds. 

Speed  regimes  where  full  annular  rub  is  either  not  in  equilibrium  or 
not  stable  may  fall  into  a  partial  rub  (if  this  mode  is  in  stable  equili¬ 
brium  at  this  rotor  speed)  or  jump  down  to  a  synchronous  whirl  without  rub. 
An  alternative  mode  of  rotor  motion  which  may  be  stable  under  these 
conditions  is  a  reverse  whirl  driven  by  the  counter-torque  produced  by  the 
frictional  forces  acting  on  the  rotor.  Some  disturbance  is  necessary  to 
perturb  the  system  into  the  reverse  whirl,  but  such  may  occur  through  a 
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transient  partial  rub  motion  following  the  deterioration  of  full  annular 
rub. 

Expanding  the  analysis  to  investigate  the  impact  of  the  stator  being 
non-rigid  results  in  extremely  detailed  equations  describing  the  motion  of 
the  system  and  the  forces  at  work  between  the  rotating  and  stationary 
element.  In  contrast  to  the  rigid-stator  system,  the  stator  now  joins  with 
the  rotor  in  responding  to  the  unbalance  driving  the  vibrations.  In 
general,  equilibrium  of  the  full  annular  rub  diminishes  in  the  vicinity  of 
the  combined  rotor-stator  natural  frequency  (which  previously  was  infinite). 
Additionally,  the  amplitude  of  the  motion  varies  greatly  now  that  the  stator 
is  not  absolutely  constraining.  If  the  stator  natural  frequency  is  greater 
than  the  rotor  natural  frequency,  then  the  recurrence  of  rub  at  speeds  above 
this  requires  some  perturbation  to  initiate  contact.  In  general,  the  finite 
dynamic  characteristics  assigned  to  the  stator  which  make  it  non-rigid 
diminish  the  reaction  forces  experienced  by  the  rotor  and  stator  (as 
compared  to  the  rigid  stator  case)  and  in  fact  may  create  a  dynamic 
relationship  less  favorable  for  full  annular  synchronous  rub. 

Much  work  remains  to  be  accomplished  in  the  study  of  rotor  rub. 
Experimental  investigation  of  the  results  presented  here  would  prove 
particularly  meaningful  in  validating  the  results  and  justifying  the 
assumptions  made.  An  analytical  investigation  of  the  equilibrium  conditions 
for  partial  rotor  rub  would  be  very  complementary  to  this  study  on  the  full 
annular  rub  problem.  Of  particular  significance  would  be  work  which 
contributes  to  the  analytical  determination  of  the  stability  of  full  annular 
rub,  partial  rub  and  reverse  whirl  for  a  rotational  system  as  a  function  of 
the  excitation  frequency. 
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APPENDIX  A 


It  is  the  intent  in  this  Appendix  to  investigate  the  dynamics  of 
several  exemplary  rotor-stator  systems  in  accordance  with  the  analysis  thus 
far  conducted  for  rotor-rub.  Unlimited  combinations  of  rotor  and  stator 
dynamic  characteristics  are  possible,  but  here  a  limited  number  are  studied 
to  demonstrate  the  impact  on  system  behavior  of  varying  key  parameters.  In 
particular,  variations  in  rotor  and  stator  natural  frequencies,  stiffnesses, 
masses  and  damping  characteristics  are  investigated.  Similarly,  the  effects 
of  varying  the  magnitude  of  the  unbalance  and  friction  are  studied.  The 
indicators  for  system  behavior  deemed  here  to  convey  the  most  information 
are: 


(i)  Rotor  Amplitude  vs.  Speed 

(ii)  Normal  Force  vs.  Speed 

(iii)  Rotor  and  Stator  Displacement  History 

Synchronous  rub  is  presumed  in  the  analysis,  and  is  indicated  by  stable 
rotor  displacement  solutions  (aside  from  the  fundamental  free-rotor 
solution)  in  conjunction  with  a  real,  positive  reaction  (or  normal)  force  at 
the  rotor-stator  interface. 

The  plot  of  rotor  amplitude  (z  /r  )  vs.  frequency  ratio  (w/w„  )  is 

c  c  Nr 

presented  for  the  two  solutions  to  Equation  (29)  as  well  as  the  free-rotor 
solution  to  Equation  (iJ).  The  free-rotor  solution  (depicted  by  the  broken- 
line)  is  the  only  stable  solution  below  the  rub  speed,  B  and  is  again 

stable  at  speeds  greater  than  B'  (described  in  section  3).  Both  equilibrium 
solutions  to  Equation  (29)  are  presented  (depicted  by  the  solid  lines 


labelled  1  and  2),  but  it  is  only  curve  1,  that  solution  initiating  at  the 
rub  speed,  which  provides  stable  equilibrium.  This  is  somewhat  intuitive  by 
analyzing  the  geometry  of  the  rotor  and  stator  equilibrium  conditions  (such 
as  Figure  23)  corresponding  to  a  particular  speed  of  interest.  Where  the 
solution  is  discontinuous,  the  recurrence  of  this  stability  is  not  assured. 

The  normal  force  is  also  plotted  for  the  stable  solution  as  a  function 
of  rotor  speed.  Regions  where  the  normal  force  is  negative  indicates  rub  is 
not  in  action.  It  is  immediately  apparent  that  no  solution  exists  for 
synchronous  rub  in  the  vicinity  of  the  combined  rotor-stator  and  the  stator 
natural  frequencies.  Stability  of  equilibrium  is  determinable  using  a 
simulation  program  provided  by  Professor  Udo  Fischer  which  traces  both  the 
rotor's  and  stator's  displacement  as  a  function  of  time.  Several  examples 
of  both  stable  and  unstable  equilibrium  are  provided  throughout  this 
appendix.  A  stable  equilibrium  is  indicated  by  the  constant  orbit  of  both 
the  rotor  and  stator.  This  solution  may  be  artificially  perturbed  (through 
initial  conditions  to  the  program),  but  still  returns  to  the  equilibrium 
condition.  Unstable  equilibrium  is  evidenced  by  varying  amplitude  of  both 
rotor  and  stator  displacement,  typically  ending  in  the  rotor  undergoing  a 
synchronous  whirl  without  rub  while  the  stator  returns  to  the  rest  position. 

Two  primary  groupings  are  analyzed:  and  <  ui^.  A  total 

of  18  cases  are  presented  with  their  key  parameters  outlined  in  Table  A1 . 
Rotor  characteristics  are  indicated  by  subscript  1  while  stator  characteris¬ 
tics  are  indicated  by  subscript  2.  One  baseline  case  for  each  of  the 
primary  groupings  is  analyzed,  with  comments  regarding  each  subsequent  case 
within  the  grouping  being  referenced  to  the  baseline  case  (unless  specifi¬ 
cally  stated  otherwise). 

These  cases,  as  described  in  the  following,  clearly  demonstrate  the 
impact  of  varying  the  system's  parameters  as  previously  discussed: 
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Case  1  This  case  provides  the  baseline  for  the  condition  .  As  an 

example  of  limited  synchronous  rub,  only  a  limited  speed  regimes  exists 
wherein  synchronous  rub  is  in  equilibrium. 

Case  2:  This  case  demonstrates  the  impact  of  an  increase  in  the  mass 
unbalance.  As  an  example  of  synchronous  rub-limited  mass  unbalance.  As  an 
example  of  synchronous  rub-limited  discontinuity,  synchronous  rub  is  in 
equilibrium  for  all  speeds  greater  than  the  rub  speed,  except  for  a  limited 
region  demarked  by  the  combined  rotor-stator  and  the  stator  natural  frequen¬ 
cies.  Additionally,  the  amplitude  of  displacement  is  increased  above  the 
baseline  while  the  reaction  force  increases  severely. 

Case  3:  A  reduction  in  friction  increases  the  speed  regime  wherein 
synchronous  rub  is  in  equilibrium.  Concurrent  with  this  is  a  large  increase 
in  the  reaction  force. 

Case  M:  Similar  to  Case  3,  a  reduction  in  stator  damping  leads  to  estab¬ 
lishing  the  equilibrium  of  synchronous  rub  over  all  speeds  outside  of  the 
combined  rotor-stator  and  stator  natural  frequencies.  As  well,  there  is  an 
increase  in  the  reaction  force. 

Case  5 ■  Reducing  the  rotor's  damping  results  in  a  large  increase  in  the 
amplitude  of  displacement  (above  the  baseline)  and  a  slight  increase  in  the 
speed  region  for  synchronous  rub.  No  recurrence  of  rub  occurs  though  above 
the  stator  natural  frequency. 

Case  6 :  An  increase  in  excitation  and  damping  on  the  rotor  and  stator 
demonstrates  how  the  damping  reduces  the  displacement  and  reaction  force  at 
low  speeds  (in  comparison  to  Case  2),  but  has  little  impact  at  high  speeds 
U  >  u^). 

Case  7:  Reducing  stator  rigidity  while  maintaining  the  same  stator  natural 
frequency  results  in  a  reduced  region  of  synchronous  rub  and  a  reduced 


reaction  force. 
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Case  8:  Converse  to  Case  7,  increasing  stator  rigidity  increases  the  region 
of  synchronous  rub  and  results  in  increased  reaction  forces. 

Case  9:  Increasing  stator  stiffness  and  so  the  stator  natural  frequency 
extends  the  speed  regime  for  synchronous  rub  while  causing  an  increase  in 
the  reaction  force  and  a  reduction  in  displacement. 

Case  10:  Increasing  stator  stiffness  and  the  excitation  demonstrates  the 
large  displacements  and  severe  reaction  forces  that  may  develop.  Once 
again,  synchronous  rub  is  re-established  above  the  stator  natural  frequency. 

Case  1 1  :  This  case  represents  the  baseline  for  the  condition  u)^s  <  A 

limited  synchronous  rub  occurs  with  a  displacement  less  than  the  radial 
clearance. 

Case  12:  Increasing  excitation  results  in  a  continuous  synchronous  rub  with 
significant  reaction  forces  developing  at  high  speeds. 

Case  13:  Reducing  friction  also  causes  a  continuous  synchronous  rub  to 
occur  similar  to  the  previous  case. 

Case  Ut:  Reducing  the  stator's  damping  has  only  a  slight  impact  on  the 
reaction  force  and  equilibrium  range  for  synchronous  rub  at  speeds  much 
greater  than  the  stator  natural  frequency. 

Case  15:  Reducing  rotor  damping  results  in  a  continuous  synchronous  rub 
similar  to  Case  12. 

Case  16:  Increasing  stator  rigidity,  while  maintaining  constant  natural 
frequency  reduces  the  speed  regime  of  synchronous  rub  in  this  case  while 
causing  the  amplitude  to  approach  the  radial  clearance. 

Case  17:  Converse  to  the  previous  case  reducing  stator  rigidity  for  a  given 
uYjs ,  causes  a  reduction  in  amplitude  with  a  slight  increase  in  the  speed 
regime  for  synchronous  rub. 
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Case  18:  Increasing  stator  mass  (and  thus  reducing  ui^)  has  little  effect 
on  the  speed  regime  or  reaction  force  of  the  synchronous  rub,  but  does  drive 
the  displacement  toward  the  radial  clearance  (i.e.,  the  rigid  stator 
condition) . 

The  stability  of  the  solutions  represented  by  displacement  curves  1  and 
2  is  analyzed  for  the  three  rub  behavior  types:  limited  synchronous  rub 
(Cases  1  and  11),  synchronous  rub-limited  discontinuity  (Case  2)  and 
continuous  synchronous  rub  (Case  12).  Within  these  examples,  several  speeds 
are  investigated  in  order  to  study  the  stability  of  each  speed  regime 
wherein  synchronous  rub  is  in  equilibrium  and  too  the  stability  of  the  two 
solutions  at  those  speeds.  These  examples  demonstrate  the  inherent 
instability  of  curve  2  solutions  as  well  as  the  variability  of  stability 
associated  with  curve  1  solutions.  In  particular,  curve  1  appears  stable 
for  both  the  limited  rub  and  continuous  rub  conditions,  but  there  is  a  speed 
range  within  the  synchronous  rub-limited  discontinuity  equilibrium  solutions 
wherein  stability  is  lost. 
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